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Abstract 

We derive explicit results from an exact expression for fermion loops with an ar- 
bitrary number of density vertices in two dimensions at zero temperature, which has 
been obtained recently by Feldman et al. |l|]. The 3-loop is an elementary function 
of the three external momenta and frequencies, and the N-loop can be expressed 
as a linear combination of 3-loops with coefficients that are rational functions of 
momenta and frequencies. We show that the divergencies of single loops for low 
energy and small momenta cancel each other when loops with permuted external 
variables are summed. The symmetrized N-loop, i.e. the connected N-point density 
correlation function of the Fermi gas, does not diverge for low energies and small 
momenta. In the dynamical limit, where momenta scale to zero at fixed finite en- 
ergy variables, the symmetrized N-loop vanishes as the (2N— 2)-th power of the scale 
parameter. 
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1. Introduction 

The discovery of high temperature superconductivity and the pecuhar phenomena ob- 
served in two-dimensional electron gases have stimulated much interest in the low energy 
physics of interacting Fermi systems with singular interactions and reduced dimensional- 
ity. In particular, the possible breakdown of Fermi liquid theory in such systems has been 
a challenging issue 0. 

In this context the properties of fermion loops (see Fig. 1) with density or current ver- 
tices play a very important role. Such loops appear as subdiagrams of Feynman diagrams 
and also as kernels in effective actions for collective degrees of freedom. Their behavior 
for small external momenta and energies determines many properties of the whole system. 
A single loop with N vertices diverges in the small energy- momentum limit for N > 2. 
In the symmetrized loop, obtained by summing all permutations of the external energy- 
momentum variables gi, . . . , gN, at least the leading contributions from single loops in the 
small energy- momentum limit are known to cancel each other 0, This systematic 
cancellation is known as loop cancellation. As a consequence, loops with N > 2 are often 
irrelevant for the low-energy behavior, and can therefore be neglected in effective actions 
or in resummations of the perturbation series via Ward identities or bosonization ^. 
This implies that collective density fluctuations and the associate response functions can 
be described by a renormalized random phase approximation, as in a Fermi liquid, even 
if interactions destroy Landau quasi-particles. 

While the 2-loop, known as polarization insertion or bubble diagram, has been com- 
puted long ago in one, two and three dimensions, no explicit formulae are available for 
higher order loops. Furthermore, the available proofs of loop cancellation show only that 
some cancellation occurs, but they do not show by how many powers the divergence is 
reduced. To be able to carry out an accurate power-counting of the infrared divergencies 
of systems with singular interactions or soft bosonic modes, a more detailed knowledge of 
the behavior of loops is required. 

Recently, Feldman, Knorrer, Sinclair and Trubowitz [|l| have obtained an exact ex- 
pression for the N-loop with density vertices in a two dimensional Fermi gas at zero 
temperature. In this article we derive explicit formulae for loops with density vertices 
in two dimensions from their expression, which are useful for the evaluation of Feynman 
diagrams containing such loops as subdiagrams. We show that loop cancellation reduces 
the degree of divergence by N— 2 powers for N > 2, i.e. the symmetrized loops, which are 
equal to the N-point density correlation functions of the Fermi gas, are generically finite 
in the small energy-momentum limit. In the dynamical limit, where the momenta scale 
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Figure 1: The N-loop with its energy-momentum labels. 



to zero at fixed finite energy variables, the symmetrized N-loop vanishes as the (2N— 2)-th 
power of the scale parameter. 

This article is organized as follows. In Sec. 2 we define single and symmetrized loops 
and review briefiy known basic properties. In Sec. 3 we summarize the main results on 
loops in two dimensions obtained by Feldman et al. which are the starting point for our 
analysis. In Sec. 4 the 3-loop is evaluated explicitly, and its asymptotic behavior for small 
energy and momentum variables is analysed. In Sec. 5 we derive a theorem clarifying the 
structure of a certain symmetrized product, which allows us to control the asymptotic low 
energy and small momentum behavior of symmetrized N-loops, to be discussed in Sec. 6. 
A conclusion follows in Sec. 7. 

2. Definitions and simple properties of loops 

The amplitude of the N-loop with density vertices, represented by the Feynman dia- 
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gram in Fig. 1, is given by 

nN(gi,...,gN) =/n(pi,...,Pn) = J Y[Go{k-pj) (1) 

at temperature zero. Here k = (/co,k), qj = {qjo,(ij) and Pj = {pjo,Pj) are energy- 
momentum vectors. The variables qj and pj are related by the linear transformation 

(lj=Pj+i-Pj, J = 1,---,N (2) 

where Pn+i = Pi- Energy and momentum conservation at all vertices yields the restriction 
gi + . . . + gN = 0. The variables gi, . . . , gN fix pi, . . . ,pn only up to a constant shift 
Pj 1-^ Pj + p. Setting pi = 0, one gets 

P2 = gi 

P3 = qi + q2 

pn = qi + q2 + ■■■ + qN-i (3) 

We use the imaginary time representation, with a non-interacting propagator 

Go{k) = ^ r (4) 

^ko - (ek - yu) 

where is the dispersion relation and yU the chemical potential of the system. The 
/cQ-integral in Eq. (0) can be easily carried out using the residue theorem; one obtains |^ 

/.to.....P„)^g/,_^,,^,^(^ n/.w| (5) 



where /ij(k) = ek-p, - ek-p^. + i{pio - Pjo). 

The 2-loop 112 (g, — g) = n(g) is known as polarization insertion or particle- hole bubble, 
and has a direct physical meaning: n(g) is the dynamical density-density correlation 
function of a non-interacting Fermi system 0. While the 2-loop has been computed 
explicitly for continuum systems in one, two and three dimensions already long ago 
only few explicit results exist for higher loops. 

The behavior of loops for small qj is particularly important, because interactions and 
bosonic propagators attached to the fermion loops in general Feynman diagrams are often 
singular for small energies and momenta H. 
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A simple formula for the static small-q limit has been derived long ago by Hertz and 
Klenin H, namely 



hm hm nN(gi, . . . ,gN) 



-"o,;- o"-^^^' • • • ' (N-l)! deN-2 



(6) 



e=ii 



where D{e) = J 5{e — e^) is the density of states. Note that the limit qj — is 
not unique; the above result is valid only in the so-called static limit where the energy 
variables go to zero first. 

In the following we will analyze the behavior of loops in the generic small-q limit, 
lim^^o nN(Agi, . . . , A^n) for arbitrary gi, . . . , ^n, with a vanishing scaling factor A applied 
to all energy and momentum variables. According to a simple power-counting estimate of 
the integral in Eq. (|I]) one would expect that an N-loop with N > 2 diverges as A^~^ for 
A — i> 0, since for qj = one integrates over an N-fold divergence on the (c? — l)-dimensional 
manifold defined by = 0, Ck = /i, which has codimension two in a (c? -|- l)-dimensional 
energy-momentum space. This estimate yields however only an upper bound for the degree 
of divergence, since the actual value of the integral may be smaller due to cancellations 
of contributions with opposite signs, as is obviously the case in the static limit (H). 

We will also analyze the so-called dynamical limit, where a vanishing scaling factor A 
is applied only to momentum variables, at fixed finite energy variables. 

Divergencies of single loops in the small-q limit cancel each other at least to some 
degree in the symmetrized loop 

^Uqi,- ■ ■ =SU^{qi,...,q^) = ^ ^Nlgpi, • • • , ^pn) (7) 

where the symmetrization operator S imposes summation over all permutations of gi, . . . , 
gN- The behavior of single loops is physically not relevant, since all physical properties 
can be expressed in terms of symmetrized loops. This important point is however ignored 
in some approximation schemes. The symmetrized N-loop is proportional to the N-point 
density correlation function of a non-interacting Fermi system with dispersion relation ek 
and chemical potential fi. In the one-dimensional Luttinger model 0, which has a linear 
dispersion relation, loops with N > 2 cancel completely, i.e. 11^ = 0, as first noticed by 
Dzyaloshinskii and Larkin ^ . In other systems the loop cancellation is not complete and, 
in general, II-^ does not vanish in the small-q limit. Only in the dynamical limit Ward 
identities associated with particle number conservation imply that 11^ goes to zero at least 



linearly in each momentum variable p|, [TO]. This is crucial in particular for the infrared 



renormalizability of the Coulomb gas |10], since the small-q singularity of the interaction 
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is not removed by screening in the dynamical limit. For the general small-q limit with 
finite ratios of momenta and energies, only the cancellation of leading singularities has 
already been established [§, H, ^ , but to our knowledge the behavior of the remainder has 
not yet been analyzed. 



3. Two dimensions and the results by Feldman et al. 

We now focus on two-dimensional Fermi systems with a quadratic dispersion relation 
ek = k^/2m, where m is the mass of the particles. The Fermi surface is then a circle 
with radius kp = \/2mfi. In the following we choose energy and momentum units such 
that kp = 1 and m = 1. It is easy to reinsert arbitrary values for kp and m by making a 
simple dimensional analysis of the expressions. 



The 2-loop in two dimensions, first computed by Stern [|TT|, is given by 



Zn 



1 

|q| 



\ 



f^Qo , |q| 



where the complex square- root is defined to have positiv real part. In the small-q limit 
one obtains 



lim Il(Xq) = — 



1 



i + |qlV%', 



(9) 

at finite 



i.e. a finite result depending on the ratio |q|/5'o- In the dynamical limit, q 
go, the 2-loop vanishes as |qp. 

Using methods from complex analysis, Feldman et al. flj have recently reduced the 
problem of computing arbitrary N-loops in two dimensions to elementary integrals. We 
now report their results, which are the starting point for our further explicit evaluations. 

If the three vectors pi, pj and p^ are not coUinear, one can define the two-dimensional 
complex vector d'^^^ as the unique solution of the equations fij(k) = fjk(X) = 0) where 



■k 



1. 2 

2 (p. 



Pj) + ^iP^O - Pjo) 



(10) 



Since fij + fjk + fki = 0, the vector d*-'^ is also a solution of the equation /fci(k) = 0. 
Hence, d*-^'^ does not depend on the order of the indices i,j, k. The real part of d*-''^ is the 
center of the circle circumscribing the triangle with vertices pj, Pj, p^. 
For the 3-loop, Feldman et al. |jl[ have obtained the expression 



hiPl,P2,P3) 



1 



27ridet(p2-pi,P3-pi^ 



dz 



'111 
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where S12 = S23 = S31 = 1, S21 = S32 = S13 = -1, and det(u, v) = Ur,Vy - UyVx for 
arbitrary two-dimensional vectors u, v. The contour-integrals are performed along the 
curves •jij = {wij{s)\0 < s < 1} where Wij{s) is the unique (generally complex) root of 
the quadratic equation 



(Pj-Pi)^ + 2 det(d-pi, Pj-Pi) z + (d-pi 
satisfying the condition 



|2 = s2 



Im '^^^"^ ~ ^^^^^^^'^^^ +dy- Piy ^ Q 
{Pjy - Piy)Wij{s) +dx - Pix 



(12) 



(13) 



with d = {dx, dy 



1123 



The general N-loop can also be written as a complex contour integral over rational 
functions [|l| . That expression implies that the N-loop can be expressed in terms of 3-loops 
as 



/n(Pi,---,Pn) = J2 

l<i<j<k<N 



n -1 



n /-(d 



ijk\ 



hiPi,Pj,Pk) 



(14) 



This reduction formula follows also directly from the identity 



Y[Go{k-pi 



1=1 



E 

l<i<j<k<N 



-I -1 



n ud''') 



Goik-p,)Goik-pj)Goik-pk) (15) 



valid for Go{k) = [ik^ — (k^/2 — jj,)] in two dimensions [Q. Note that fi^id 



ijk\ 



ijk\ 



fkv{d 



ijk\ 



A remarkable explicit result has been obtained by Feldman et al. [jT| in the static case 
{Qio = for all i). In that case they show that In vanishes identically for N > 2 when 
all disks with radius kp around the points pi, . . . , pn have at least one point in common! 
This implies that nN(q'i, . . . , Qn) vanishes identically for = and small q^, if N > 2, 
which is compatible with Eq. (^) since D{e) is constant in two dimensions. 

In the following we will derive more explicit formulae for 3-loops, and will analyze the 
behavior of IIn and 11^ in various important limits. To this end, it is useful to determine 
the point d^^^ explicitly as a function of Pi,Pj,Pk- Writing the inhomogeneous linear 
system of equations for d^^^ in matrix form and inverting the matrix, one obtains 



^(Pfc - P'i) + i{Pko - Pio) 



(Pi 



+ j ^ k (16) 



det(pj-pi,pA,.-pi 

Here p-*- = {px,Py)^ = {—Py,Px) denotes a rotation by an angle tt/2 and j ^ k means 
exchange of indices j and k. Note that |det(pj— pj, Pfe— Pi)| is twice the area of the triangle 
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with vertices Pi,Pj,Pk- If Pi,Pj,Pk are collinear, this area vanishes and the point d^^^ 
tends to infinity. Using (]T6|), the coefficient fiu{d'^^'^) appearing in the reduction formula 
(|l^) can be expressed exphcitly as 

1 



+ < ^{pI - p1) + KPko - Pio) 



det{pj-pi,p^-pi 
det(pj— Pi,pfc-pi 



+ j 



(17) 



We conclude this section by defining a complex function of Pi,Pj,Pk that will be useful in 
the following evaluations: 



Inserting d*-''^, one obtains 
\Pj-Pi\ 



det(p,-d^^■^Pi-d^J■'= 



Pi 



iPj-Pil 



det d^^'^-p,, 



Pj-Pi 
\Pj-Pi\ 



•^ijk 
Uijk 



2det(pj-pi,pfc-pi 



(Pi-Pfc) ■ (Pfc-Pi) 



Pj-Pi 



det(pj-pi, pfc-pi) |pj-pi 



• [iPk-Pi)iPjo-Pko) - iPj-Pk)iPko-Pio)] (19) 



The modulus of Xijk is twice the distance between the point Re d*-^*^ and the straight line 
connecting pi and pj. Note the useful identity 

2 



\2 -2 



^ijk 



1, I , . PjO-PiO 

2 Pi- Pi 



(20) 



4. The 3-loop 

In this section we derive an explicit formula for the 3-loop and analyze its behavior in 
various important limits. 

We first derive an explicit expression for the functions Wij{s), i,j G {1,2,3}, i ^ j. 
Equation (|T2|) can also be written as 

[|p,-p,|z + z,,f + (d-p,)2 - 4 = (21) 

where Zij = Xij + yij = Zijk with k G {1, 2, 3}, k ^ Making the ansatz 

(S) = (22) 

\pj- Pil 
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with a complex function %(s) = Xij{s) +iyij{s), and using relation (^^, one obtains the 
equation 



+ 



1, I , . PjO-PiO 

o Pi-Pi +«i r 

2 Pj-Pi 



(23) 



Splitting real and imaginary parts, one obtains two real equations for Xij{s) and yij{s) 



Xijis) - Vijis) 
Xij{s)yij{s) 
and the condition (|13l) implies 



IPj-Pil + 



\pj-Pi\^ 



aij[s) 



1, 



iPjo-Pio)xijis) - -\pj-pi\ yij{s) > 



(24) 



(25) 



The two equations (^) define intersecting hyperbolas in the {x, |/)-plane. Note that 
the second equation does not depend on the variable s. For any s there is exactly one 
intersection point satisfying the condition {^^), given by 



Xij[S) 



sgn{pjo-pio)-^ \j \l[aij{s)Y + {Pjo-Piof + aij{s) 



(26) 



All the square-roots in ( pUf ) have real positive arguments. 

It is now obvious that the curve = {wij{s)\Q < s < 1} is connected and sweeps out 
an angle smaller than vr with respect to the origin in the complex plane. Hence 



(27) 



/ ^ = ln 


'Wijil) 







where In is the main-branch of the complex logarithm, which is real and continuous on 
the positive real axis and has a branch cut on the negative real axis. Hence, 



hiPl,P2,P3) 



2711 det(p2-pi,P3-p 



E 111 



%■(!) 



1; i,j=i 



:28) 



We have thus obtained an explicit formula for the 3-loop as a function of Pi,P2,P3, which is 
useful in particular for the evaluation of Feynman diagrams with 3-loops as subdiagrams. 
Note that the result J3 = in the static limit (for momenta such that unit disks around 
them have at least one point in common) is not obvious from Eq. (pH]), and is obtained 
only as a consequence of cancellations in the sum over 
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We now analyze the behavior of single 3-loops and the symmetrized 3-loop in the 
small-q limit defined in Sec. 2. With the choice pi = 0, the substitution qi i-^ Xqi implies 
Pi I— > Xpi. The function yij is invariant under this substitution, while Xij becomes a 



homogeneous function of A of order one, i.e. Xij 
and yij{s) reduce to 



Xxij. For A — > 0, the functions Xij{s) 



Xij[S) 



sg: 



n{Pjo-Pio)\J s'^ + (Pjo-Pio)VIPj 



Pi 



A 



I Pi -Pi 



2 a/I + s2|pj-p,|2/(pjo-pio)^ 



(29) 



Hence Xij{s) is of order one for A — > while yij{s) is of order A. To avoid confusion, we 
emphasize that we consider the small-q limit for finite fixed ratios between the various mo- 
menta and frequencies. The following small-q expansion is valid only if energy /momentum 
ratios f-'""^''! etc. do not vanish. Hence, one cannot recover the static limit from the 
asymptotic expressions. The expansion does hold, however, in the dynamical limit, where 
energy /momentum ratios tend to infinity. 

Using the power-series expansion ln(l — z) = — J2'^=i ^-2", we can expand 



In 



In [xij{s) 



Wij\ - 



00 1 



n=l 



n 



(30) 



s=0 



The n-th term in the sum is of order A"". To compute J3, we have to sum over all pairs 
with i,j G {1,2,3} and i 7^ j. Since 



X 



) Vji ~ ~yij and Xji{s) — —Xij{s) , yji{s) — yij{s) 



(31) 



the first term in (|30|) cancels when the contributions from and 7^, are subtracted (see 
Eq. (^), i.e. 



dz 

7«j Z 



dz 

7ji Z 



00 -I 



n=l 



n 



1 


•^ij ~l~ Wij ( "5 ) 


n 


•^ij ^Vij ( ) 


1 




Xjj(s) iyij _ 




X ( S ) iyij _ 





(32) 



s=0 



The leading term for small A is the one with n = 1. Neglecting the other terms, one 
obtains after some elementary algebra (see Appendix A), 

1 1 



hiPl,P2,P3) A" 



X 



27ridet(p2-pi,P3-pi) l + (Imd)2 
^ det(Imd, pj -pi) 

(ij)=(l,2),(2,3),(3,l) + \Pj-Pi?/{PjO-PiQy 



(33) 



with corrections of order one. Hence, a single 3-loop generally diverges as A ^ for A ^ 0, 
just as expected from simple power-counting. 
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We now determine the symmetrized 3-loop for small energies and momenta. There are 
only two non-equivalent permutations of qi,q2,<l3 (see Fig. 2), corresponding to a change 
of sign of pi,p2,P3. The symmetrized loop Ilf can therefore be expressed as 



nf(gi,g2,g3) = ^ [h{PuP2,P3) + h{-Pl,-P2,-P3)] = ^^(Pl,P2,P3) 



(34) 



It is easy to see that 



Xi 



X. 



Vij ^ -Vij and Xij{s) ^ -Xij{s) , yij{s) ^ yij{s) for pi h-> -p^. (35) 



Hence, the first order terms cancel when summing permuted loops, and the leading con- 
tributions come from n = 2 in the expansion (p2D, i.e. 



nf(gi,g2,g3) 



where 



A-ni det(p2-pi,P3-pi) 



E 

(i,i)=(l,2),(2,3),(3,l) 



KM 



•^ij i ^Uij ("5) 



After some simple algebra (see Appendix A) one obtains 

1 1 



i=0 

(36) 
(37) 



n3(gi,g2,g3) 



E 



X 



[s" + (Imd) 



7rdet(p2-pi,P3-pi) [s2+(Imd)2]2 (• 2) (2,3),(3,1) 

|Pi-P*l/2 



Xij 



1 + s'^\pj-pi\'^/{pjo-p^o) 



- \Pjo-Pio\ 



\ 



jPjo - Pio) 

Ipj-Pip 



(3^ 



s=0 



with corrections of order A. The symmetrized loop is thus finite (of order one) for A — 0, 
i.e. the divergence present in single loops is cancelled when summing permutations. Note 
also that Hf (gi, q2, gs) is obviously real in the small-q limit. 

Next we analyze the behavior of the symmetrized 3-loop for one vanishing external 
momentum variable, say qi, while the other momenta and all energy variables remain 
finite. Noting that det(p2 — pi, P3 — pi) is of order qi in that limit, it is easy to see that 
X^{s) = 0(|qi|) for all Hence, we can use the expansion (^) once again. For the 



symmetrized loop the first order terms cancel, and the leading contribution is given by 
(^), which is of order |qi|. 

We finally consider the dynamical limit, where all momenta scale to zero (with a 
scaling factor A) at fixed finite energy variables. In that limit X^(s) is of order A^, so 
that the leading contribution to the symmetrized loop is again given by (pG]). Solving for 
Zij{s) for small A, one can easily show that the difference X^{1) — X^(0) is only of order 
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A*^. Since det(p2— pi, Ps— Pi) is obviously of order in the dynamical limit, one obtains 
the result 

nf((gio,Aqi),(g2o,Aq2),(g3o,Aq3)) =0(A^) for A^O (39) 

5. Permutations of external energy-momentum variables 

Consider a term in the reduction formula (|14D with fixed values i,j, k. Let {%, qi, . . . , 
Qm} C. {qi, . . . ,qiy} be the subset of those external energy-momentum variables that are 
attached to the loop between the lines labelled by pj and pk, as shown in Fig. 3 (the case 
of variables attached between pi and pj or p^ and Pi is completely analogous). Define 
Po = Pj, Pi = Pj+i, P2 = Pj+2, ...,Pm= Pk-l, Pm+1 = Pk- 

Our aim is to show that the sum 

M ^ 

P fi=i J fj. 

over all permutations P of qivi is of order one in the small-q limit, and vanishes 

with a power 2M in the dynamical limit. Here 

= /fl+,(d^^') = (P? - P.) ■ d^'' + IIpI - (pJ)'] + ^iP^o - P%) (41) 

is constructed with the energy-momentum variable p^ on the /i-th line resulting from the 
permutation P. Note that for a fixed permutation P, the product n^li(l//^) diverges 
as A~*^ for A ^ 0. Only the sum over all permutations turns out to be finite in that limit, 
due to systematic cancellations. 

As a first step, we sum over permutations leading to sequences 

gi, gs, • • • , go, ga+i, • • • , gM, (42) 

where gi, . . . , g^ follow the order of their indices, while go is placed between ga and ga+i- 
We denote these particular permutations by (a), where a G {0, 1, ... , M}. Obviously 

= I = ^M+i - ^0 for /i < a ^^^^ 
^ \ p^ for /i > a 

Define 

= fl%M'') = (pi^^ - p^) ■ d'^' + llp" - (p1"^)'] + 'iP^o - & 

/; for /i<a ^^^^ 
for fi> a 
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i i 



12 
i 



1m-i 1m 

i i 



k-p. k-pj k-p^ k-pj 



Figure 3: Part of an N-loop from pj to pk- 



i.e. ffj, is constructed with and with p'^ = Vfi + % ~ Qo- For the sum over all 
permutations (a) we have derived the following 



Lemma 1 : 



E n = 7 E qo ■ (qi + • • • + qa) n 7H 



(45) 



) /x=l J a=l 

Proof: Using p'^ = p^+i — pi + Pj and the identity /o = fij{&^'') = 0, it is easy to derive 
the relation 

/i + /; = /m+i+^m+i (46) 



where 



This yields 



ga+i = qo ■ (pq+i-pi) = qo ■ (qi + • • • + q«) 



M 

./'l fa. M=l fi 



1^1 



/g+l + 9a+l 
fl 



M 

n 



(47) 



(45 



We now prove the identity (]45|) by starting from the right hand side. Relation ( ^8|) yields 



M 

E 



1 /l ;t^L /m 



M M 



M 



M 



n 



{a) 



En7T^+TE(/:-/«+i)n 



=1 ^t=l ^1 «=l 



1 f^" 



(a) 



(49) 



The second term reduces to n^=i(l//Ai)) due to a cancellation of all terms but the first 
one in the a-sum (note that /a/+i = /ifc(d*-^'^) = 0). This completes the proof of Lemma 
1. For energy- momentum variables attached between pi and pj or between p^ and p^, it 
can be proved in the same way. 

For the following formal manipulations it will be helpful to represent products of 



to the permutation P, i.e. 



inverse /-factors X^^=\{^U^) by the ordered sequence of momentum labels corresponding 



M 1 



= {P0,Pl,P2,...,PM) 



(50) 



14 



With this notation, Lemma 1 can also be written as 

M M 

J2{l,2,...,a,Q,a + l,...,M) = j (qo ■ Qp) (1, 2, . . . , a, 0, a + 1, . . . , M) (51) 

P<q: 

In particular, for M = 1 one obtains 

(0,1) + (1,0) = ^^(1,0) (52) 

and for M = 2, 

(0, 1, 2) + (1, 0, 2) + (1, 2, 0) = ^ {(qo • qi) (1, 0, 2) + [qo • (qi + qs)] (1, 2, 0)} (53) 

The sum over permutations (a) has reduced the degree of divergence for A ^ by 
one power. We now consider permutations {a,P) with {3 <a leading to sequences 

^2, ^3, ■ ■ ■ , QP+U ■■■,qa, %, Qa+l, ■ ■ ■ , QM (54) 

i.e. q2, ■ ■ ■ Am follow the order of their indices while qi is placed at an arbitrary position 
before go in the sequence. The corresponding energy-momentum variables on fermion 
lines become 

Pi = -%-qi for ii<(3 

= p^+i - go for P < ii<a (55) 

^ pfj, for > a 

and the associated product of inverse f-factors is 

fi-^ =(2,3,...,/3,l,/3 + l,...,o;,0,a+l,...,M) (56) 

The sum over all permutations {a, (3) with fixed a and (3 < a, can be rewritten using 
Lemma 2: 

En7(i^ = 77Eqi-(q2 + ... + q/3+i)n7(i^ + tzIIt;; IT t ^^^^ 

The second term on the right hand side can be written symbolically as 

1 a— 1 \ ^ \ 1 

-JJ^-fT, n :^ = T7(2,3,...,a,[0,l],a+l,...,M) (58) 

where [0, 1] stands for the sum of two energy-momentum variables, qo + Qi, i-e. the product 
of inverse f-factors above corresponds to the sequence q2,Q3, ■ ■ ■ ,Qa,Qo + Qi, Qa+i, ■ ■ ■ Am- 
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Proof: The proof of Lemma 2 is almost identical to the one for Lemma 1. A difference 
arises only in the last step, where now the "boundary term" (^) contributes, because 
/q, 7^ replaces the factor /m+i = (see Eq. (^91) and below). 

For example, for M = 2, a = 1, Lemma 2 yields the trivial identity 

(1,0,2) = 1 ([0,1], 2) (59) 
Jl 

and for M = 2, a = 2, 

(1, 2, 0) + (2, 1, 0) = ^ (qi ■ qs) (2, 1,0) + ^ (2, [0, 1]) (60) 
Jl Jl 

An example for M = 3 can be found in Appendix B. 

Lemma 1 and Lemma 2 have been formulated and derived for go and qi "running" 
from left to right, with a specific fixed order of the other variables. Analogous identities 
hold of course for any permutation of running and fixed variables, e.g. for ^3 running with 
some fixed order of all other variables. In general, if qp is the running variable, the explicit 
prefactor 1/fi in Lemma 1 and l//( in Lemma 2 has to be replaced by a factor 1/hp, 
where 

hp = (Pi + qp - Pi) ■ d'^'' + 2 ~ (Pj' + ^P^^^ + " (^JO + ^"o)] (61) 

i.e. an /-function constructed with pi and pj + qp. Note that /i = ho and f[ = hi by 
definition. 

We are now ready to prove the main result of this section: 

Theorem 1: The sum SMiqo, qi, ■ ■ ■ , qn) = J2pllff=i Tp can be written as a sum over 
fractions with numerators 

(q^i ■ q^O (q<x2 ■ q^) ... (q^,, ■ q^;,) (62) 
where 7^ cr^, and denominators 

M-l 

ho hi ... hM X Yl_ /-factors (63) 
1 

where the f -factors are constructed with pi andpj + some variables q. In each numerator 
each momentum variable from {qo, . . . , (\m} appears at least once. 

Proof: The following algorithm reduces the sum Yl,pT{^u=ijp to the form described in 
the theorem. We proceed by induction with respect to the number of energy-momentum 
variables M + 1. We will prove that the algorithm works for the case M = 1 directly, and 
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for general M under the assumption that it works for less than M + 1 variables. 

Step 1 : Write down all permutations of (0, 1, ... , M), and apply Lemma 1 to groups with 
a fixed order of 1, ... , M, and running over all possible positions. This yields 

En^ = Enr'E(---,A---,o,...) (64) 

where, for fixed p, all products (. . . , p, . . . , 0, . . .) with p to the left of contribute exactly 
once. For M = 1 one gets simply 

Siiqo, qi) = (0, 1) + (1, 0) = (1, 0) = (65) 

tlQ tlQ til 

1. e. for this case Theorem 1 is already proven. 

Step 2: The coefficients (. . . , p, . . . , 0, . . .) of can be further reduced by applying 
Lemma 2 to groups with fixed positions of all variables except p and p running from the 
first place to the place before and nearest to 0. This yields a sum 

E ^^E(---.P'.---.A---.0,...) + 7^E(---J0,p],...) (66) 
p'=l '^p '^p 

where the coefficient (. . . , p', . . . , p, . . . , 0, . . .) is a sum over all permutations with p' to 
the left of p and p to the left of 0. The "boundary term" 5] (. . . , [0, p], . . .) is a sum over 
all permutations of products of M — 1 /-factors with M energy-momentum variables, i.e. 
all variables with p 7^ 0, p and % -\- qp. According to our induction hypothesis we can 
apply Theorem 1 to this term, which, with the prefactor obtained in Step 1 and Step 

2, leads to a sum over fractions with the required structure. 

Step 3: To reduce the first sum of terms obtained in Step 2, we apply Lemma 2 to the 
coefficients (. . . , p', . ■ ■ , P, ■ ■ ■ , 0, . . .), this time with p' running from the first place to 
p. This yields sums 

E ^^E(---'P^---'P''---'A---'0,...)+-^i:(...,[p,p'],...,o,...) 

p"¥=o,p,p' 

(67) 

Again, the boundary term leads back to a case with one variable less, which is reducible 
to the desired form by virtue of the induction hypothesis. 

Steps 4-M: The first sum of terms from Step 3 can be reduced by another application 
of Lemma 2, and so on. The boundary terms always lead back to a case with a reduced 
number of energy-momentum variables. The algorithm terminates after the M-th step. 
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where coefficients {p^^~^\ p^^~'^\ . . . , p", p', p, 0) with a fixed order appear, where p(^~^^ 
is the only variable that has not "run" (in the sense of Lemma 2). These coefficients 
correspond to inverse products of M /-factors, the first of which is hp(M-i). It is clear 
that each momentum variable appears at least once in each product of M scalar products 
generated by the algorithm. This completes the proof of the theorem. 

For M = 2, for example, the algorithm yields 



S2{qo, qi, q-2) = (qo ■ qi) (qi • q2^ 



ho hi /i2 



1 1 



L^o,i hi; 



cychc permutations (68) 



where hp^p' is an /-function constructed with pj and pj + qp + qp'- In Appendix B the 
algorithm is presented at work for the case M — 

The following important corollaries follow directly from Theorem 1. 

Corollary 1 : 

SM{qo, qi, ■ ■ ■ , qu) is finite (of order one) and real in the small-q limit, where all momen- 
tum and energy variables scale to zero as A. 

Corollary 2: 

Sniqa, qi-i ■ ■ ■ Am) is of order for qp — > if g^o and all with n ^ p remain finite. 
Corollary 3: 

In the dynamical limit, where all momenta scale to zero as A at fixed finite energy vari- 
ables, the sum 5'm(?0) qi, ■ ■ ■ , ?m) vanishes as X"^^ . 

Proof: Theorem 1 tells us that 5'm(?o, ■ ■ ■ , ?m) is a sum of fractions with M scalar 
products of momenta in the numerator and 2M /-factors in the denominator. The /- 
factors vanish linearly in the small-q limit (where momenta and energies tend both to 
zero), and are purely imaginary to leading order in A. This yields Corollary 1. The 
/-factors remain finite if momenta vanish at finite energies. This yields the Corollaries 
2 and 3, where for the former the observation that in each numerator each momentum 
variable appears at least once is crucial. 

Actually CorroUary 2 holds already for a sum over all positions of qp with fixed order 
of the other variables (instead of summing all permutations), and follows directly from 
Lemma 1. 

In the following sections we will use these results to control the behavior of sym- 
metrized N-loops in the above limits. 



6. The N-loop 
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Any N-loop can be computed by inserting the explicit expression for the 3-loop given 
in Sec. 4 into the reduction formula (p^, with /^^(d*-''^) from (]T7|). 

In the small-q limit, where all momenta and energy variables scale to zero as A, the 
unsymmetrized N-loop diverges as A^~^, since J3 diverges as and the product of N— 3 
factors /ij/(d*-'^) vanishes as A^~^. We will now show that the symmetrized N-loop (|^) 
remains finite in the small-q limit. 



Symmetrizing the reduction formula (|Tj) for the N-loop, one can write 



HnI^i, • • • ,9n) = 5 Sij SjkSkih{Pi,Pj,Pk) (69) 

l<i<j <k<N 

Here S'jj, Sjk and S^i are symmetrized products of inverse /-factors, e.g. 

I 1 for k = i + 1 

with Sm as defined in Sec. 5. Note that symmetrizing once or twice, or first symmetrizing 
partially (with respect to a subset of variables) and then completely (by applying 5), or 
vice versa, always yields the same result. In particular, 

n^(gi, . . . ,gN) = ^ [n^(gi, . . . ,gN) + (Pi, • • • ,Pn) (-Pi, • • • , -Pn)] (71) 

Since fiu{d^^'') ^ fl{^'^^) and thus Sm ^ ^Jf for (pi, . . . ,pn) ^ (-pi, • • • , -Pn), one 
can also write 

n^(gi, . . . ,gN) =S Yl [Re{SijSjkSki)li{puPj,Pk)+i'^'oci{SijSjkSki)I^{Pi,Pj,Pk) 

l<i<j<k<N 

(72) 

where li{pi,Pj,Pk) = \h{phPj,Pk) - hi-Pi, -Pj, -Pk)- Note that unsymmetrized loops 
are complex for N > 2, while symmetrized loops are real, since /n(— Pi, . . . , — Pn) = 
^n(Pi>--->Pn)- 

In the small-q limit, Re(S'jj Sjk Ski) is of order one, while Im(S'jj Sjk Ski) is of order A 
(see Theorem 1 and Corollary 1 in Sec. 5). In Sec. 4 we have shown that is of order 
one and real in that limit, while diverges as A~^ with a purely imaginary coefficient 
(see Eq. (133|)). Hence H§(gi, . . . , ^n) is finite in the small-q limit, i.e. 

H^(Agi,...,AgN) = 0(l) for A^O (73) 

We now analyse the dynamical limit, where all momenta scale to zero as A at finite 
energy variables. In that limit Re/j;^(d*-'^) is of order A^, while Im/jj,(d*-'^) is of order 
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one (see, for example, Eq. (|T7p. Theorem 1 then imphes that Re{SijSjkSki) vanishes as 
_)^2(N-3)^ and Im(SijSjkSki) as A^'-^"^-'. The symmetrized 3-loop has been shown to vanish 
as A^ in the dynamical limit in Sec. 4. Since X^{s) is of order A^ and - Xj(0) 

of order A^ (see Sec. 4), it is obvious from (|32|) that vanishes as A^ in the dynamical 
limit. Inserting these asymptotic power-laws in (|72|), one finds 

n^((gio,Aqi),...,(gNo,AqN) = 0(A2N-2) for A^O (74) 

We finally discuss the behavior of loops for a single vanishing momentum, say qi, 
while energy variables and other momenta remain finite. Note that this limit cannot be 
taken for 2-loops, since momentum conservation imposes q2 = — qi in that case, i.e. q2 
cannot remain finite if qi tends to zero. Unsymmetrized loops with N > 3 remain finite in 
the small qi-limit, since /-factors and 3-loops remain finite if the other variables remain 
finite. We have already seen in Sec. 4 that the symmetrized 3-loop vanishes linearly for 
qi —>■ 0. This behavior holds also for the symmetrized N-loop, i.e. 

ng(gi,...,gN) = 0(|qi|) for qi ^ (75) 

This latter result can be established quite easily in any dimension by using the techniques 
reviewed in Refs. 0, ^, for example via Ward identities 0, |10 . 



The above results for the symmetrized loops have been checked numerically for N 
3, 4, . . . , 7. No cancellations beyond the order established here have been observed. 



7. Conclusion 

In summary, we have derived explicit formulae for fermion loops in two dimensions 
from an exact expression obtained recently by Feldman et al. . The 3-loop is an elemen- 
tary function of momenta and frequencies, and the N-loop can be expressed as a linear 
combination of 3-loops with coefficients that are rational functions of momentum and fre- 
quency variables. These formulae are very useful for the evaluation of Feynman diagrams 
for two-dimensional interacting Fermi systems where such loops appear as subdiagrams. 

We have also analyzed to what extent divergencies in the low energy and small mo- 
mentum limit cancel in the symmetrized loop, defined as a sum over permutations of 
momentum and frequency variables. A single loop was seen to diverge with a power N — 2 
as expected from a simple scaling analysis. We have proved a reduction of the degree of 
divergence in the symmetrized N-loop by N — 2 powers, such that symmetrized loops and 
thus all N-point density correlation functions of the two-dimensional Fermi gas are gener- 
ally finite in the low energy and small momentum limit. This drastic cancellation justifies 
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the frequent neglect of loops with more than two insertions assumed in Ward identity or 
bosonization methods [0, Q for systems where interactions with small momentum trans- 
fers dominate the low-energy physics. In the dynamical limit, where momenta scale to 
zero at finite energy variables, we have shown that the symmetrized N-loop vanishes as 
the (2N — 2)-th power of the scale parameter. 
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Appendix A: 3-loop for small energies and momenta 



In this appendix we present the algebra leading to the asymptotic results for the 3-loop 
in Sec. 4. 



The first order term in (|3^) is 



2iyij{s 



s=l 



s=0 



2i 



[xij{s) +iyij]yij{s) 



(76) 



s=0 



The denominator can be simplified to 



(Imd) 



Imd 



Pj-Pi 
\Pj-Pi 



+ (Imd)^ 



(77) 



where in the second step we have used the fact that d is a zero of the function /jj(k) 
defined in (0). Relation ( p4|) yields the identity 



Xij{s)yij{s) = Q 

(i,j)={l,2),(2,3),{3,l) 



and yijyij{0) = |det(Imd, p^ -p^) yields 



y^jy^j{0) = 

(ij)=(l,2),(2,3),(3,l) 



(7J 



(79) 



Inserting the above relations in the first order term in (|32D one immediately obtains the 
result 



The second order term in (|32D is 

s=l 



2i Xij yij (s) 



2i 



Xij yij{s)[xjAs) - yfj + 2ixij{s) y^j] 



s=0 



[s2 + (Imd)2]2 



s=l 



^0) 
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where the last step follows from ([77|). Using Eq. (0), i.e. Xij{s) yij{s) = —^{pjo —Pio), 
and the identity 

Vij (pjo - pio) = (81) 



(ij)=(l,2),(2,3),(3,l) 



one obtains the result 



Appendix B: Permutations for M = 3 

In this Appendix we carry out the permutation algorithm constructed in Sec. 5 for 
the case M = 3. 

Step 1: Write down all permutations of (0, 1,2,3), i.e. 24 terms, and apply Lemma 1 to 
groups with fixed order of 1, 2, 3. This yields 



go • qi 

ho 



:i, 0, 2, 3) + (1, 2, 0, 3) + (2, 1, 0, 3) + (1, 2, 3, 0) + (2, 1, 3, 0) + (2, 3, 1, 0) + 2 ^ 3 



f82) 



and analogous terms proportional to ^7^^ and . 

Step 2: Applying Lemma 2, the above coefficients of can be rewritten as 

1 



(1,0,2,3) 
(1,2, 0,3) + (2, 1,0, 3) 
1,2,3,0) + (2, 1,3,0) + (2, 3, 1,0) 



hi 

qi ■q2 



([0,1], 2, 3) 



1 



(2, 1,0, 3) + -(2, [0,1], 3) 
III III 

■ (2, 1, 3, 0) + • + (2, 3, 1, 0) 



hi 

+ 7^ (2, 3, [0,1]) 
hi 



hi 



^3) 



and analogously for permutations 2 3. Applying Lemma 1 to the "boundary terms'' 
yields 

1 

h'l 



([0,1], 2, 3) + (2, [0,1], 3) + (2, 3, [0,1]) 

(qo + qi) ■ q2 



(2, [0,1], 3) 



(qo + qi) ■ (q2 + qs) 



hiho^i hiho^i 
and analogously for permutations 2^3. Hence, the coefficient of becomes 



(2, 3, [0,1]) (84) 



qi • q2 

hi 

(qo + qi) ■ q2 



hihoA 



(2, 1, 0, 3) + (2, 1, 3, 0) + (2, 3, 1, 0) + (3, 2, 1, 0) 
(2, [0,1], 3) + (2, 3, [0,1]) + (3, 2, [0,1]) 



+ 2^3 



^5) 
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The coefficients of ^"^'^^ and ^"^^^^ can be rewritten similarly by applying Lemma 2 to 
groups with 2 and 3 running from the left end to the position of with the order of all 
other variables fixed. 

Step 3: We finally show, as an example, how another scalar product of momenta can be 
extracted from the coefficient (2, 1, 0, 3) + (2, 1, 3, 0) + (2, 3, 1, 0) + (3, 2, 1, 0) of ^ 
obtained in the second step above. Lemma 2 yields 

(2,1,0,3) = -^([1,2], 0,3) 

'^2 



(2,1,3,0) 
(2, 3, 1,0) + (3, 2, 1,0) 



1 



h2 

q2 ■ qa 



([1,2], 3,0) 



(3,2,1,0) + -(3,[1,2],0) 
h2 h2 



Applying Lemma 2 to the boundary terms in (|86|) yields 



(86) 



([1,2], 0,3) 



([1,2],3,0) + (3,[1,2],0) 



(qi + q2) • qs 



,1,2], 3) 



(3, [1,2],0) + 



1 



For the boundary terms above. Lemma 1 yields finally 

1 r,.. . ,^ . (qo + qi + q2) ■ qa 



(3, [0,1, 2]) (87) 



([0,1, 2], 3) + (3, [0,1,2]) 



^2^1,2^0,1,2^3 



Thus only terms which contain products of three scalar products of momenta are left over. 
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